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Azimuthal angular correlations of charged hadrons with respect to the axis of a reconstructed 

(trigger) jet in Au+Au and p+p collisions at ^snn = 200 GcV in STAR are presented. The trigger 
jet population in Au+Au collisions is biased towards jets that have not interacted with the medium, 
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allowing easier matching of jet energies between Au+Au and p+p collisions while enhancing medium 
effects on the recoil jet. The associated hadron yield of the recoil jet is significantly suppressed at high 
transverse momentum {p?^^°'^) and enhanced at low in Au+Au collisions, which is indicative 

of medium-induced parton energy loss in ultrarelativistic heavy-ion collisions. 

PACS numbers: 25.75.-q, 25.75.Bh, 25.75. Gz, 12.38.Mh, 21.65. Qr 



High-energy collisions of heavy nuclei at the Rela- 
tivistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Laboratory produce an energy density for which a 
strongly-coupled medium of deconfined quarks and glu- 
ons is expected to form [11-0]. The properties of this 
medium can be studied using partons with large trans- 
verse momenta (pt) resulting from hard scatterings in 
the initial stages of the collision. The scattered partons 
recoil and fragment into back-to-back clusters of hadrons, 
known as jets. 

Jets in p+p collisions are well-described by pQCD the- 
ory [5] and can be used as a reference for studies of 
medium-induced jet modification. Comparisons are also 
frequently made between jets in Au-|-Au and rf+Au in 
order to separate initial and final state effects [B] E]- 
Suppression of high-px hadrons in single-particle mea- 
surements [8HTT] and of particle yields on the recoil side 
("awayside") of high-px triggered "dihadron" correla- 
tions P^^HTB] has been observed in more-central Au-f Au 
collisions relative to p+p and d+Au collisions. This sup- 
pression of jet fragments is often attributed to partonic 
energy loss due to interactions with the medium [19j . 
However, these single- and two-particle methods do not 
allow direct access to the original parton kinematics. Fur- 
thermore, due to the steeply falling parton energy spec- 
trum, the trigger hadron population in dihadron correla- 
tions is dominated by fragments from lower-energy par- 
tons. 

Due to recent advances in jet-finding techniques in 
heavy-ion collisions [50] it is now possible to study trig- 
gered correlations with respect to the axis of a recon- 
structed jet, instead of using the dihadron correlation 
technique in which a high-px hadron is used as a proxy 
for the jet axis. Jet reconstruction allows more direct ac- 
cess to the original parton energy and makes it possible to 
select a sample of higher-energy partons, thus increasing 
the kinematic reach of these correlation measurements. 

Conceptually, models of partonic energy loss fall into 
two general categories: (1) radiative/coUisional energy 
loss models [121 f^T^ES] , in which partons lose energy and 
are scattered as they traverse the medium (radiative en- 
ergy loss is expected to dominate for light partons), and 
(2) calculations within a gauge/gravity duality frame- 
work, specifically in a strongly-coupled Af — A supersym- 
metric Yang-Mills theory [2BH55], in which light partons 
either emerge from the medium nearly unmodified or de- 
posit all their energy into the medium. In this analysis, 
azimuthal angular correlations of mid-rapidity charged 
hadrons are studied with respect to a reconstructed mid- 



rapidity (trigger) jet. In the radiative energy loss model 
the correlations should exhibit a difference in the widths 
and/or yields of the jet peaks in Au-|-Au compared to 
p+p, while in AdS/CFT models little or no jet shape 
modification is expected (although jet suppression would 
be observed in the awayside associated particle yields 
and in single-particle inclusive measurements). The ef- 
fects of medium-induced partonic energy loss, or "jet- 
quenching," can be studied by comparing the shapes and 
associated hadron yields of jets in Au-|-Au with those in 
p+p collisions. 

The data used in this analysis were collected by the 
STAR detector at RHIC for p+p and Au-|-Au colli- 
sions at ^snn = 200 GeV in 2006 and 2007, respec- 
tively. Charged tracks are reconstructed in the Time 
Projection Chamber (TPC) [29] and the transverse en- 
ergy (Et) of neutral hadrons is measured in the Barrel 
Electromagnetic Calorimeter (BEMC) towers (with size 
X Ary = 0.05 x 0.05) l30j. A 100% hadronic correc- 
tion scheme (in which the total momentum of all charged 
tracks which point towards a BEMC tower is subtracted 
from the energy of that tower) is applied to account for 
charged hadron energy deposition in the BEMC. BEMC 
towers associated with electrons are rejected. 

Events are selected by an online high tower (HT) trig- 
ger, which imposed thresholds on £'x of about 5.4 GeV 
to be deposited in at least one BEMC tower. A com- 
mon offline software HT threshold of > 6 GeV is 
imposed (after hadronic correction). In Au-|-Au only the 
20% most central events are analyzed, where event cen- 
trality is determined by the uncorrected charged particle 
multiplicity in the TPC within pseudorapidity [Tyj < 0.5. 
Events are required to have a primary vertex position 
along the beam axis within 25 cm of the center of the 
TPC. Tracks are required to have pr > 0.2 GeV/c, at 
least 20 points measured in the TPC (out of a maximum 
of 45) , a distance of closest approach to the collision ver- 
tex of less than 1 cm, and |?7| < 1. Events containing 
tracks with px > 30 GeV/c are not considered in this 
analysis because of poor momentum resolution. Particle 
distributions are corrected for single particle tracking ef- 
ficiency and for detector pair acceptance by event mixing 
(in relative azimuthal angle A(f> only). 

Jets are reconstructed from charged tracks in the TPC 
and neutral towers in the BEMC using the anti-Zc-r al- 
gorithm from the FastJet package PJ] with a resolu- 
tion parameter R — 0.4. The reconstructed jet axis is 
required to be within \ri\ < 1 — R. Only tracks with 
Px > 2 GeV/c and towers with i?x > 2 GeV are used 
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FIG. 1. (Color online.) (a) j^^^''""" spectra of HT trigger jets 
in p+p and Au+Au, and of p+p HT trigger jets embedded in 
Au+Au MB events, (b) The ratio of dN/dj^^' to 
dAf/dpi^'"'^'"^'*"''''^". Uncertainties due to the relative tracking 
efficiency, relative tower energy, and AE shift are shown as 
shaded bands. Note that pi^''^'^'^ is calculated only from tracks 
and towers with pt > 2 GeV/c. 



FIG. 2. (Color online.) Jet-hadron correlations are shown 
after background subtraction for 10 < pi^^''^°'^ < 15 GeV/c 
and for two ranges in pSf""'^: (a) 0.5 < pT"" < 1 GeV/c and 
(b) 4 < Pt'°'' < 6 GeV/c. The data points from Au+Au and 
p+p collisions are shown with Gaussian fits to the jet peaks 
and systematic uncertainty bands due to: tracking efficiency, 
the shape of the combinatoric background, and the trigger jet 
energy scale defined by the AE and SDaa shifts. 



in the jet reconstruction in order to control the effects of 
background fluctuations. Furthermore, the reconstructed 
jet must include as one of its constituents a BEMC tower 
that fired the HT trigger. While in most jet reconstruc- 
tion analyses it is necessary to subtract an average back- 
ground energy from the reconstructed jet |32j . the 
2 GeV cut on tracks and towers makes this subtraction 
unnecessary here. Instead, background fluctuations are 
handled using an embedding procedure described later. 

Jet-hadron correlations are defined as distributions in 
Acj) = 0jot — (/'assoc, where (fijct denotes the azimuthal 
angle of the axis of a reconstructed (trigger) jet and the 
associated particles are all charged hadrons, measured 
as TPC tracks, in the event. To obtain the associated 
particle yields (Y) and widths (cr) of the jet peaks, the 
raw correlation functions are flt with the functional form: 

+ B (l + 2i;r°'=4" cos(2A0) + 2v^''''°''vf' cos(3A0)) , 



which includes two Gaussians representing the trig- 
ger/nearside (NS) and awayside (AS) jet peaks, and a 
background term modulated by V2^^°'^V2^ and v^^^'^'^ir'^^ . 

The elliptic anisotropy of the background is assumed to 
factorize into the product of the single-particle anisotropy 
of the associated particles due to elliptic flow (w|®''°'^) and 
the correlation of the jet axis with the 2"'^-harmonic event 
plane [53] ('^2'''). The possibility that there is a corre- 
lation between the jet axis and the 3'^'^-harmonic event 
plane |34j (which can give rise to a Wg"* component and 
thus a non-zero vf^^°'^v'^^ term) is also taken into account. 
The values of V2^^°'^v''2*' and v^^^°'^u'^^ are discussed later. 

The Gaussian yields of the jet peaks, Y, are integrated 
over a given bin in the transverse momentum of the as- 
sociated hadrons (p^^""^), and the reconstructed jet px 
(p:!^*'™), as well as over the A77 acceptance. 

The effects of medium-induced modification can be 
quantified by the widths of the jet peaks, cr, as well as 
Daa and EDaa, defined in Eqns. ^ and ([3| respec- 
tively. Daa measures the transverse-momentum differ- 
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FIG. 3. (Color online.) The Gaussian widths of the awayside 
jet peaks {a as) in Au+Au (solid symbols) and p+p (open 
symbols) are shown for two ranges of p!^^'""^: 10 — 15 GeV/c 
(red circles) and 20 — 40 GeV/c (black squares). The results 
for 15 — 20 GeV/c (not shown) are similar. The boundaries of 
the p^^°'^ bins are shown along the upper axis. YaJEM-DE 
model calculations (solid and dashed lines) are from [35] . 

ence between Au+Au and p+p (in a given p^^°^ bin with 
mean {p^^°^)): 

DAAiPTn =yAu+Au{pTn ■ {pt°1au+au (2) 

V (n^Bsoc\ / assocN 
-fp+plPx ) \Pt /p+p- 

YiDaa nieasures the energy balance over the entire p^^°'^ 
range: 

SZ?AA= DAAipT"")- (3) 

If jets in Au+Au and p+p have identical fragmentation 
patterns, then Daa = for all Deviations from 

Daa = are indicative of jet modification. 

In order to make meaningful quantitative comparisons 
between jets in Au+Au and p+p, it is necessary to com- 
pare jets with similar energies in the two collision sys- 
tems. While the reconstructed jet pr is not directly 
related to the original parton energy (especially in this 
analysis because p^r^'''^"'^ is calculated only from tracks and 
towers with > 2 GeV/c), jets in Au+Au with a given 
y^t.rec, Au+Au rnatched to similar p+p jets using the 
following procedure: The effect of the background associ- 
ated with heavy-ion collisions on the trigger jet energy is 
assessed through embedding p+p HT events in Au+Au 
minimum bias (MB) events (with the same centrality and 
high-multiplicity bias as the Au+Au HT events). Under 
the assumption that Au+Au HT trigger jets are simi- 
lar to p+p HT trigger jets in a Au+Au collision back- 
ground, the correspondence between the p+p jet energy 
(pj^*'™'P+'') and the Au+Au jet energy (^^^^''^'^'^^p+p ^"^^ ~ 
y^t,rec,Au+Au-j determined through this embed- 

ding. Figure [l] compares the pi^^''^°'^'P~^P '^"^^ spectrum 
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FIG. 4. (Color online.) The awayside momentum difference 
Daa is shown for two ranges of 10 — 15 GeV/c (red 

circles) and 20 — 40 GeV/c (black squares). The results for 
15 — 20 GeV/c (not shown) are similar. The boundaries of 
the p^^°'^ bins are shown along the upper axis. YaJEM-DE 
model calculations (solid lines) are from [35j . 

to the pj^*'''<='='Au-i-Au gpg(,|;].^jjj measured in Au+Au HT 
events. For a given range in j)J^''''°'^'P+p the correspond- 
ing distribution is obtained. When compar- 
ing Au+Au jets to equivalent p+p jets in this analysis, 
the Au+Au signal is weighted according to this distribu- 
tion. This procedure largely accounts for the effects of 
background fluctuations in Au+Au events; the possibil- 
ity of additional discrepancies between the reconstructed 
jet energies in Au+Au and p+p will be included within 
systematic uncertainties described below. 

The performance of the TPC and BEMC can vary in 
different collision systems and over time. These varia- 
tions are accounted for in the relative tracking efficiency 
between Au+Au and p+p (90% ±7% for px > 2 GeV/c), 
the relative tower efficiency (98% ±2%), and the relative 
tower energy scale (100% ± 2%). These variations in de- 
tector performance were included, and their systematic 
uncertainties were assessed, in the p+p HT (g) Au+Au 
MB embedding. The effects of the relative tracking effi- 
ciency uncertainty and the tower energy scale uncertainty 
on the spectrum are shown in Fig. [ijb). The 

embedding also accounted for jet V2 and its associated 
uncertainty (discussed later) by weighting the distribu- 
tion of the p+p HT jets with respect to the event planes 
of the Au+Au MB events; different hadronic correction 
schemes were also investigated. The effects of the tower 
efficiency and jet V2 on the jet energy scale are found to 
be negligible, as is the effect of the hadronic correction 
scheme on the final results. 

In order to analyze the jet correlation signal in Au+Au 
collisions it is necessary to subtract the large combina- 
toric background in heavy-ion collisions. The background 
levels are estimated by fitting the functional form in ([l]) 
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iet.rcc 

Pt 


T.DAA 


Detector 


V2 and V3 


Jet Energy Scale 


(GeV/c) 


(GeV/c) 


Uncertainties 


Uncertainties 


Uncertainties 


10-15 


-0.6 ±0.2 


+0.2 
-0.2 


+3.7 
-0.5 


+2.3 
-0.0 


15-20 


-1.8±0.3 


+0.3 
-0.3 


+ 1.0 

-0.0 


+1.9 
-0.0 


20-40 


-1.0 ±0.8 


+0.1 
-0.8 


+ 1.2 

-0.1 


+0.3 
-0.0 



TABLE I. Awayside EDaa values are shown with their statistical uncertainties as well as uncertainties due to detector effects, 
the shape of the combinatoric background, and the trigger jet energy scale (given by the AE and EDaa shifts). 



to the raw Acf) distributions in Au±Au and p+p, with the 
flow terms constrained to zero for the latter. The shape of 
the Au±Au background is not well-constrained because 
v^2^ and U3°* have not yet been measured experimentally. 
Therefore the uncertainties are investigated using two di- 
ametrically opposed assumptions. To assess the effect of 
the uncertainty in the shape of the background, the as- 
sumption is made that Au±Au HT trigger jets undergo 
no medium modification. Then to assess the effect of 
the uncertainty in the jet energy scale, the assumption is 
made that Au±Au HT trigger jets are maximally modi- 
fied, which is the case when Uj"* = ^3^* = 0. 

First, it is assumed that the Au+Au HT trigger jets 
undergo no modification and are equivalent to p+p HT 
trigger jets (at all p^^°^). It is expected, based on jet 
quenching model studies [35] , that the combination of the 
constituent cut and the HT trigger requirement biases 
the jet population towards unmodified jets. When fitting 
the A(j) distributions with the functional form in ([T]) the 
nearside yields and widths in Au±Au are fixed to the 
values measured in p+p, V2^^°'^v^2^ is fixed to a reason- 
able mean value and v^^^°'^v^;^^ is left as a free param- 
eter. The mean V2^^°'^ is estimated to be the average 
of z;2{FTPC}(p^"''°=) and t;2{4}(p^""°'=), while 4°' is esti- 
mated to be z;2{FTPC}(6 GeV/c), where wslFTPC} and 
U2{4} are parameterized as functions of from STAR 
minimum bias data in |37) . Here, W2{FTPC} is estimated 
with respect to the event plane determined in the For- 
ward Time Projection Chambers (2.4 < |?7| < 4.2) [35] 
and W2{4} is determined using the 4-particle cumulant 
method [32] ■ The vf^^°'^v'^^ values that result from the 
fits are reasonable when compared to measurements from 
PHENIX [3D]. The systematic uncertainties are deter- 
mined by fixing V2^'^°'^vi^^ to maximum and minimum 
values while letting v^^'^°'^v^^^ float to force the Au±Au 
nearside yields to match p+p. The lower and upper 
limits on V2^^°'^ are estimated to be V2{4:}{p^^°'^) and 
W2{FTPC}(p^""°=). The lower and upper bounds on wf* 
are conservatively estimated to be 70% and 130% of 
?;2{FTPC}(6 GeV/c), respectively. Additionally, it is ob- 
served in Fig.[l]that the shape of the jet energy spectrum 
in Au± Au does not quite match the spectrum of p+p HT 
jets embedded in Au±Au MB events, even after account- 
ing for detector effects. The spectrum shape mismatch is 
attributed to an incomplete understanding of the back- 



ground. The difference in spectra shapes is encompassed 
by a AE = +1 GcV/c shift in the Au+Au trigger jet pt, 
as shown in Fig. [Tj 

The second assumption is that the Au+Au HT trigger 
jets are maximally modified compared to p+p HT trigger 
jets. The background conditions that allow maximum in- 
creases in the nearside widths and yields are V2'^^°'^v''2^ = 
and v^^^°'^vf*' — 0. Under this assumption, SDaa — 
when the parent parton energies are correctly matched, 

,1 1 ict,rcc, Au+Au / ict,rGC,r>+p i ict,rcc ■ 

even though pij, 7^ Ft because plj, is 

calculated only from charged tracks and neutral towers 
with Pt > 2 GeV/c. The shift in the Au+Au trigger 
jet energy necessary to force SDaa to zero defines an- 
other systematic uncertainty estimate. Figure [2] shows 
background-subtracted A(f) distributions with systematic 
uncertainties due to the two opposing assumptions de- 
scribed above. 

The nearside jet is expected to have a surface bias [41] 
which makes it more likely that the recoil parton will 
travel a significant distance through the medium, there- 
fore enhancing partonic energy loss effects and quenching 
on the awayside. The awayside Gaussian widths (ctas), 
-Daa, and SDaa are reported here for different ranges 
in The awayside widths, shown in Fig.jsj at high 

^assoc g^j^g ^-j^g same in p+p and Au+Au on average, in- 
dicating that jets containing high-px fragments are not 
largely deflected by the presence of the medium. The 
widths at low are indicative of broadening. How- 

ever, as the low-p|?''°'^ widths are anticorrelated with the 
magnitude of v^^^°'^v'^^, measurements of w^f* are nec- 
essary before quantitative conclusions are drawn. The 
awayside I?aa, shown in Fig. |4j exhibits suppression of 
high-p^'^''°'^ hadrons and enhancement of low-p^^°'^ jet 
fragments in Au+Au, indicating that the jets in Au+Au 
are signiflcantly softer than those in p+p collisions. The 
SDaa values, shown in Table |l| indicate that the high- 
^assoc suppression is balanced in large part by the low- 
^assoc enhancement. 

Theoretical calculations from YaJEM-DE [35]) a 
Monte Carlo model of in-medium shower evolution, are 
also shown for ctas and Da A in Figs, [s] and |4] [3S]. This 
model incorporates both radiative and elastic energy loss, 
and describes many high-px observables from RHIC. Af- 
ter the intrinsic transverse momentum imbalance, k-r, of 
the initial hard scattering was tuned to provide the best 
fit to the p+p baseline yields {Yas,p+p), this model of 
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jet-medium interactions largely reproduces several of the 
quantitative and qualitative features observed in data. 
At high p^^°'^ the Au+Au and p+p widths match and 
the jet yields are suppressed, while the missing energy 
appears as an enhancement and broadening of the soft 
jet fragments. 

To conclude, jet-hadron correlations, which are az- 
imuthal angular correlations triggered by reconstructed 
jets, are used for the first time to investigate the modifica- 
tion of the properties of partons that traverse the medium 
created in heavy- ion collisions. The trigger/nearside jet 
sample is highly biased towards jets that have not in- 
teracted with the medium, which may enhance the ef- 
fects of jet-quenching on the recoil/awayside jet. While 
the widths of the awayside jet peaks are suggestive of 
medium-induced broadening, they are highly dependent 
on the shape of the subtracted background; measure- 
ments of ^2°* a-nd Ug"* are needed before additional quan- 
titative conclusions about broadening can be made. It 
is observed that the suppression of the high-px asso- 
ciated particle yield is in large part balanced by low- 
possoc enhancement. The redistribution of energy from 
high-px fragments to low-px fragments that remain cor- 
related with the jet axis is quantitatively consistent with 
the YaJEM-DE model of in-medium shower modifica- 
tion, and in general, qualitatively consistent with radia- 
tive/collisional energy loss models. 
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